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Abstract. We analyze the coherence properties of neutron wave packets, after 
they have interacted with a phase shifter undergoing different kinds of statistical 
fluctuations. We give a quantitative (and operational) definition of decoherence and 
compare it to the standard deviation of the distribution of the phase shifts. We find 
that in some cases the neutron ensemble is more coherent, even though it has interacted 
with a wider (i.e. more disordered) distribution of shifts. This feature is independent of 
the particular definition of decoherence: this is shown by proposing and discussing an 
alternative definition, based on the Wigner function, that displays a similar behavior. 
We briefly discuss the notion of entropy of the shifts and find that, in general, it does 
not correspond to that of decoherence of the neutron. 

PACS numbers: 03.65.Bz, 03.75.Be, 03.75.Dg 
1. Introduction 

Decoherence is an interesting phenomenon, related to the long-standing issue of 
irreversibility. Nowadays, it discloses challenging perspectives in the light of new 
technologies and related physical applications. There is a widespread consensus ||, ^ 
about the meaning of decoherence, viewed as the loss of quantum mechanical coherence 
of a physical system in interaction with other systems ("environment"). However, a 
quantitative definition of decoherence is subtle and involves conceptual pitfalls 0]. In 
addition, it always depends on the experimental configuration. An interesting quantity 
in this context is the square of the density matrix |0. Apart from lacking idempotency 
for mixed states, this quantity enjoys other interesting features p, but also yields results 
which are at variance with naive expectations based on entropy [Q. 

In this article we will consider two different definitions of decoherence: the first is 
operational and stems from an analysis of the visibility in quantum (as well as classical, 
as we will see) interference experiments. We stress that these experiments are routinely 
performed in neutron optics [0, ^, |^ . The second definition is based on the idempotency 
defect of the density matrix and is, in this sense, less operational. 

In both cases, decoherence displays an "anomalous" behavior, both as a function 
of the features of the fluctuations and the incoming state. Some concrete examples will 




Figure 1. Scheme of a Mach-Zehnder interferometer. 



be considered and discussed. Our analysis will focus on neutron optics and hinge on 
an approach based on the analysis of statistical fluctuations [llO| , |ll|. However, since 



our results are just a consequence of the wave nature of neutrons and their coherence 
properties, we expect that the same general conclusions be valid for other quantum (and 
classical) waves. 



2. Fluctuations in neutron optics 

Let us start our analysis by considering a neutron beam that crosses a Mach-Zehnder 
interferometer (MZI), as schematically shown in Fig. ||. A phase shifter A is placed in 
the lower arm of the interferometer and lipm) is the initial wave packet. 

We neglect wave-packet dispersion effects, so that the outgoing states in the 
ordinary and extraordinary channels read 



1 + l^in), 

1 - ei^^l 



IV^in), (1) 

respectively. We focus on the ordinary channel, the analysis for the extraordinary one 
being identical. Define the operator 

0(A) = ^ [l + ei^^] , (2) 

that accounts for the state evolution in the ordinary channel, and consider the output 
density matrix 

Po = |^o)(^o| = 0(A)|^i,)(V'in|0(A)t = 6(A)pi„6(A)t, (3) 

where is the density matrix of the incoming state. The trace of po yields the relative 
frequency of neutrons in the ordinary channel. 

Suppose now that the phase shift A fluctuates according to a probability law 
w{A — Aq), Aq being the average phase (operationally defined as the phase that is 
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measured — or inferred — in an interferometric experiment). Therefore one has 

J dA w{A) = 1, J dA w{A)A = 0. (4) 

The trace of the average density matrix is 

Tr = Tr J dA w{A - Ao)0(A)pi,0(A)t = Tr (^pi,6(A)tC)(A)) , (5) 

where the bar denotes the average over the distribution w{A — Aq). One obtains, after 
some algebra, 



pA 



0(A)tO(A) = -(^l + cos^j. 

Consider now the Fourier transform of the probabihty density of the fluctuations 
Q{p) = J w{A)er^P'^ 

pA 



(6) 



pA 

dA w{A) cos + i 



dA w{A) sin ■ 



h 



= C{p)+zS{p), 

where C and 5* are respectively the real and the imaginary part of Q 

C{p) = Re n{p), 
S{p) = Im Q{p). 

In Eq. we can write 



pA f pi\ 
cos — — = / (iA w(A — Aq) cos — — 
n J n 



pA 



cos 



CiP) 



sm 



(7) 

(8) 
(9) 

S{p). (10) 



h h 

In this paper, for simplicity, we will always consider symmetric distribution functions, 
that is w{A) = w{—A). Therefore 

s{p) = 0, c{p) = m (11) 

and d^) becomes 



0(A)tO(A) 



1 + Q{p) cos 



pAp 

h 



(12) 



We notice, incidentally, that the same results are obtained with a different setup §|: 
consider a polarized neutron that interacts with a magnetic field perpendicular to its 
spin. Due to the longitudinal Stern-Gerlach effect [0, its wave packet is split into two 
components that travel with different speeds and are therefore separated in space. After 
a projection onto the initial spin state, the resulting final state is slightly different from 
that considered in the preceding equations: we need to replace l^po) (and analogously 
Hje)) in i) with 

iV'o) — |^o) = 0'(A)|V^i„), (13) 

where 



6'(A) = i 



(14) 
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and A is in this case the spatial separation between the two wave packets corresponding 
to the two spin components. By averaging over A it is easy to show that one obtains 



again (|T2]). 

By plugging the average operator (|T2D into (|) one finally gets 



Tr po=2 



fi(]5) cos 



pAo 



(15) 



where (■■•) = Tr [pin ■ ■ ■] denotes the expectation value over the initial state pin- On the 
other hand, the momentum distribution is easily shown to be 



Po{p) = {p\pd\p) = Tr(|p)(p| p^) = -PM 



1 + Q{p) cos 



pAo 

h 



(16) 



where 



PM = (pIPinb). (17) 

We now introduce the visibility of the interference pattern (in the ordinary channel) 

Po{p)mAX - Po{p) mh 



Vip) 



18) 



Po{p)max + Po{p) mil 

where Po{p)max [Po{p) min] is the maximum [minimum] value assumed by Poip) when 
Ao varies. By the very definition (|^, one can verify that < V(p) < 1. Notice that, 
according to this definition, the visibility is a function of momentum p and yields a 
measure of the fringe visibility of a postselected beam of momentum p as a function of 
the phase shift Aq 0, H . Equivalently, it is a measure of the "local" spectral visibility, 
under the assumption of a slowly varying wave envelope, and so it corresponds to (the 
absolute value of) the amplitude of the cosine function in (p!6|). By using (|^) and (p!8|), 
one infers that the visibility is the modulus of the Fourier transform of the distribution 
of the shifts A and is therefore a quantity that is closely related to the physical features 
of the phase shifter. In this way we can easily relate the visibility of the interference 
pattern (and, as we will see below, the decoherence) to the "environmental" fiuctuations. 
Note that a completely equivalent definition of the spectral visibility (|TBp, which is 
nevertheless more symmetric and makes use also of the extraordinary channel, reads 

pAo 



V{p) 



\Po{p) ~ Pe{p)\ 
max ^ , , 
Ao Po{p) 



Pe{p) 



max 

Ao 



f2(p) cos ■ 



h 



where the momentum distribution of the extraordinary channel is given by 



Pe{p) = \pUp) 



fl{p) COS 



h 



(19) 



(20) 



whence Po{p) + Pe{p) = Pm{p)- The spectral visibility in the form (^9|) leads to 
a straightforward generalization which is at the basis of an operational definition of 
decoherence. 
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3. An operational definition of decoherence 



Let us endeavor to give a quantitative definition of decofierence based on tfie definition of 
visibility given in the previous section. We start from the relative frequency of particles 
detected in the ordinary and extraordinary channels 



Tr Po = ^ 



Trp^ = - 



Their difference is 

and one can define a generalized visibility 
V = max |Aro(Ao) - AfEi^o) 



1 + (^^{p) COS 
1 — ( ^l{p) COS 



pAo' 

n 
h 



^e{Ao) = ( n{p) cos 



pAo' 



(21) 



(22) 



max 

Ao 



r2(p) cos 



pAo' 

h 



max 

Ao 



dp Pi^{p)fl{p) cos 



(23) 



It is apparent that (p3D is the straightforward generalization of the spectral visibility 
(pjp, because obviously A/q+A/'e' = 1. It represents a global feature of the outgoing state, 
in contrast with the local character of (0). Notice, however, that when Pin(p') = S{p'—p) 
(incoming monochromatic beam of momentum p), the generalized visibility (^) reduces 
to the standard "local" visibility (|T8|) 



V 



max 

Ao 



max 

Ao 



dp' 6{p' — p)Q{p') cos ■ 



h 



fl{p) cos 



V{p). 



(24) 



This is a consistency check, because a spectral postselection is equivalent to injecting 
an incoming monochromatic beam. 
In general one gets 

pAo 



V < max 

~ Ao 



dpPUp)Mp)\ 



cos ■ 



dp PMV{p) = {V{p)) . (25) 



The generalized visibility yields the maximum "distance" between the intensities Mo and 
Me and is bounded by the "local" visibility averaged over the momentum distribution 
of the incoming state. 

Notice that -Pin(p) is a nonnegative quantity, while fi(p), being a Fourier transform, 
is not. For this reason, in general, the max^o does not simply enter the integral in (|23|) , 
so that (^) is a strict inequality. However, in the particular case ^^(p) > 0, V saturates 
its upper bound, the equality sign holds in (^5]), and Eq. ( p3|) simplifies into 

V = jdp PUpMp)- i^ip) > 0) (26) 

As is often to be expected, the most interesting cases are those situations in which Q{p) 
is not always positive, giving rise to "anomalous" situations. 
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Figure 2. Relative frequencies of neutrons detected in the ordinary Afo and 
extraordinary A/fi channel versus fcoAo (fco ~ po/h), for an incoming Gaussian wave 
packet (^8|) with fco*^ = 12 and a fluctuation-free phase shifter. The two intensities 
differ in phase by tt and their sum is 1. The generalized visibility (03) is 1. 



In order to understand the physical meaning of the generahzed visibihty, it is useful 
to look at the example of a fluctuation-free phase shifter, w{A) = 5(A), for which (|^) 
yields fi(p) = 1, so that the generalized visibility ( p3D becomes 



V 



max 

Ao 



dp Pin{p) COS ■ 



dp Pinip) 



for any incoming distribution Pin- This result follows also directly from (l26|). 
instance, for an incoming Gaussian wave packet 



exp 



26\ 



Po) 



(27) 
For 

(28) 



one gets the interference patterns Ao and A/g shown in Figure 0, where it is apparent 
that V = 1. 

If, on the other hand, the phase shifter fluctuates, the amplitude of the envelope 
function decreases and V < 1. We therefore give an operational definition of 
decoherence, by defining a decoherence parameter. 

pAo' 



e = l-V 



1 — max 

Ao 



max 

Ao 



r2(]3) cos ■ 



h 



dp Pin(p)fi(p) COS 



pAc 



h 



(29) 



Notice that, by Eq. (p7|), e = for a fluctuation-free phase shifter (quantum coherence 
perfectly preserved), while e ^ 1 when the magnitude of the fluctuations increases, 
fl{p) —>■ and the envelope function in Figure || squeezes away all oscillations, eventually 
yielding A/o(Ao) = Me{Ao), independently of Aq. Observe also that V and e are 
independent of the coherence of the initial state (namely, they do not depend on the 
off-diagonal terms of the density matrix). On the other hand, they strongly depend on 
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the momentum distribution of the initial state (pT]). In this sense they measure the loss 
of quantum coherence caused by a given physical setup, independently of the coherence 
of the incoming state. 

It is important to stress that the above definition of decoherence is operational. 
One first measures the relative frequencies of neutrons detected in the ordinary and 
extraordinary channels as a function of Aq, both being measurable quantities. Then 
one evaluates (|23|) and computes e. 



4. Some examples 



The decoherence parameter (|29D depends on the product of the momentum distribution 
of the incoming beam times the spectrum of the phase-shifter fluctuations, -Pin(p) x ^(p)- 
These two ingredients affect e at the same level. Therefore, their role can be 
interchanged: by maintaining their product unaltered, there exist "dual" situations 
that give exactly the same decoherence parameter with very different kinds of statistical 
fluctuations and incoming states. 

By keeping the above remark in mind, it is interesting to look at some particular 
cases that can be treated analytically. Let the phases be distributed according to a 
Gaussian law with standard deviation a 



w{A - Ao 



1 



exp 



(A-Ao 



so that = exp (^—p^a'^/2h^^ and the decoherence parameter reads 



(30) 



1 — max 

Ao 



dp P\n{p) exp 
I- j dpPin{p)exp 



p^a"^' 



cos 



V h 



p^a'^^ 



2r 



For the Gaussian wave packet (|28|) one gets 



1 



52 



with fco = Po/^- This is exact and is shown in Figure |^. At fixed 5 the decoherence 
parameter (^2|) increases with o", although the details of its behavior are strongly 
dependent on the spatial width of the packet 5. This behavior is in agreement with 
expectation: decoherence e increases with the magnitude a of the fluctuations. 
For a monochromatic beam [Pin(p') = 5{p — p')] 

£fe = l-e- — , (33) 

with k = p/h. This is shown in Figure ^(a) and can be obtained from ( p2D in the 5 ^ 00 
limit. Notice that high momenta are more fragile against fluctuations |]TT|]. Moreover, 
when the distribution of the shifts is Gaussian, Sk and equivalently V(p) are monotonic 
functions: they both depend "smoothly" on a. 



6' 



21,2' 



(31) 



(32) 
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Figure 3. Decoherence parameter e ( p^ ) versus width S of the Gaussian wave packet 
and standard deviation a of the fluctuating shifts (fco = po/ti). 



Let now the phase shifts be distributed according to the law 



u;(A- Ao) = - 



1 



TT 



2(t2 - (A - Ac 



(34) 



,1/2 



for |A — AqI < V2a and otherwise, with standard deviation (/ AP'w{/\)d/S.) ' = a. 
From an experimental perspective this is more convenient and easier to reproduce than 
the Gaussian distribution (pOl): indeed, (0) follows from a phase A(t) = Ao + -\/2crsint, 
where t ("time") is a parameter, uniformly distributed between and 2it, namely 
w{A) = Jq^ dt 6{A — y/2a sint) /2t[. (One can require y/2a < Aq, in order that A{t) 
be positive — and the term v^asint be regarded as a "small" fluctuation around the 
average value. However, strictly speaking, this is not necessary from a mathematical 
point of view.) From (|3^) and (^ one gets 



V2a clA 

dt ( 

— exp 

^7r/2 TT 



e fi 



A2 



I — - — smr 



(35) 



where Jq is the Bessel function of order zero. The decoherence parameter ([29| ) reads 



max 

Ao 



dp P\n{p)Jo 



( V2pa\ fpAo 



cos 



V h 



and for a monochromatic beam one obtains {k = p/Ti) 



1 — max 

Ao 



/ ^/2pa\ f pA, 



cos 



h 



1 - \Jo{V2ka) 



(36) 



(37) 



This function is shown in Figure H(b): observe that decoherence is not a monotonic 
function of the noise a in (041). 
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(a) (b) 

Figure 4. (a) Decoherence parameter Sk (^) versus fccr, for a monochromatic beam 
interacting with a shifter fluctuating according to (pO|); (b) Decoherence parameter 
£fc (|3^) versus ka, for a monochromatic beam interacting with a shifter fluctuating 
according to (^. 



A comparison between Figures Ka) and |5(b) is interesting. In both cases one 
observes fragility at high momenta p = hk. However, the behavior of decoherence 
in Figure ^(b) is somewhat anomalous and against naive expectation. For a given k, 
there are situations where decoherence e decreases by increasing the strength of the 
fluctuations a. Note also that we are considering incoming monochromatic beams, 
whence, according to (p4D and (p9D , = 1 — V{hk) and the decoherence parameter is 
strictly related to the standard visibility of the interference pattern. Therefore, in the 
anomalous regions, one observes an increase in visibility by increasing the fluctuations 
of the phase shifter, a phenomenon somewhat similar to stochastic resonance . This 
is true not only for monochromatic beams, but also for narrow distributions (packets) 
in momentum space. 

These anomalous results are not entirely surprising, if one compares them to other 
known results in classical optics. We will therefore recall in the next section some 
notions related to the visibility of a classical interference experiment: the visibility can be 
expressed as the Fourier transform of the spectral distribution of a quasi-monochromatic 
light source and it displays some "anomalies" even in cases that are different from our 
"Gaussian" example (|28|) . 

5. A classical analogy 

The phenomena analyzed in the previous sections have an interesting classical 
counterpart that is worth looking at in some detail. In this section we will examine 
the behavior of the visibility in a two-beam interference experiment, in relation to the 
spectral density distribution of the source. We follow Born and Wolf |T^. Suppose to 
have two beams whose optical difference is AS and whose wave number is k = 2'k/X. 
Their phase difference reads 

6{k,AS) = kAS, (38) 
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and, assuming that they have the same intensity \{k)dk in the range [k,k + dk], the 
intensity at the screen due to the elementary wave number range dk reads 

i{k, AS)dk = 2\{k)[l + cos{kAS)]dk. (39) 

Observe that the different spectral components add incoherently, so 

I{AS) = 2 j dk \{k)[l + cos{kAS)] (40) 

is the intensity at the screen as a function of AS, due to both interfering beams. The 
quantity k is to be compared to the phase A in Sec. ||. 

In some cases one deals with light sources that emit with characteristic spectral 
lines. If we consider only one of these spectral lines, \{k) is different from zero only in a 
very small range of k about some mean value /cq. Putting 

j{k) = 2\{ko + k), (41) 

the intensity at the screen ( ^OD becomes 

I{AS)= J dk j{k){l + cos[{ko + k)AS]} 

= N[1 + C{AS) cos{koAS) - S{AS) sin(A;o A5)] , (42) 

where is a normalization factor, defined as the sum of both the (equal) intensities 
of the beams, and C and S are the average value on the spectral distribution j{k) of 
cos(/cA(S) and sm{kAS) respectively 



N = Jdkj{k), (43) 
C{AS) = ^ J dk j{k)cos{kAS), (44) 

S{AS) = j^jdk j{k) sm{kAS), (45) 

i.e. C and 5* are respectively the real and the imaginary part of the Fourier transform 
Vl{AS) of 3{k)/N 

C{AS) = Re Q{AS), S{AS) = Im Q{AS), (46) 

Q{AS) = Jdk^^ e'^^^. (47) 

From (|i2|), the intensity at the screen can be written as 

/ = [1 + \Vl{AS) I cos(A;oA5 + ^{AS))] , (48) 

where taxiip{AS) = S{AS)/C{AS). A comparison with Eq. (plGj) shows that AS plays 
the same role of p. 

Because j{k) is very peaked about k = 0, variations of C and S can be considered 
negligible compared with cos(A;oA5) and sin(/coAiS) in Eq. (^2|); analogously for ip in 
(1^). Consequently, under the assumption of slowly varying envelope, one can define a 
"local" visibility [l^, given by 



V(A5) = 'itti"" I HA 9!"'° = • (^5' 
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J V J ' V 




V=|sin(£A^A5)|_ 



|iA/tA5| 




(c) 



J =Jo^ 
V~e~ 
with 



(4^)"|cos(|a5)|, 



a^k 2 



AS 




i;ioi,|i 



(d) 



V ~ 1 £.'(45") Vs + 4 cos(^ as), 
with 



aAk 



Figure 5. Visibility versus optical path difference AS, for a quasi-monochromatic 
light source with spectral distribution j{x) (pl|). Different spectral distribution shapes 
are considered: square- like (a), Gaussian (b), double Gaussian with peaks at the same 
level (c) and at different levels (d). In (b), (c) and (d) Ak = 2-\/ln2/Q;. (Reproduced 
with permission from Ref. p^.) 



expressed as a function of the optical path difference AS. The visibihty is therefore the 
amphtude of the cosine function in Eq. (^S)). Observe that, whenever j{k) is an even 
spectral distribution, 

V(A5) = |C(A5)| = |1](A5)| (50) 

and it is possible to determine (apart from the sign) the Fourier transform of j{k)/N 
from the visibility. 

Equations ([l8|) and (^91) are easily compared. The distribution of the phase shifts in 
the quantum case is replaced by the spectral distribution of the incoherent light source 
in the classical case. Indeed, f2(AiS), C{AS) and S{AS) in (|46|)-(^7D correspond to 
fl{p), C{p) and S{j)) in i.e. the Fourier transform of j{k)/N corresponds to that 
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of w{A) (notice that w{A) is normalized to unity). 

The visibihty curves ( |50D are shown in Fig. ^ for different shapes of the spectral 
distribution j{k). As one can see, they show different behavior. In Fig. ^(a) a square-like 
spectral distribution gives rise to a visibility function \ sin y/y\, in Fig. ^(b) a Gaussian 
spectral distribution produces a Gaussian visibility function, in Fig. ||(c) and (d) two 
"double Gaussian" distributions (with the peaks that have or do not have the same level, 
respectively) yield more complicate visibility functions. Only in the case (b), i.e. with a 
Gaussian spectral distribution, the visibility is a monotonic function of the optical path 
difference AS. In such a case, the naive expectation is confirmed that, by increasing 
the optical path difference, the visibility decreases. This is not true in cases (a), (c) and 
(d), where the visibility is not a decreasing function for every range of AS, but there 
are regions on the screen where, by increasing the optical path difference, the two-beam 
interference visibility increases. 

Similar results can be obtained if one considers two-beam interference with extended 
monochromatic light sources. In such a case, the source is treated as a collection of 
monochromatic point-like sources that add incoherently and, instead of j{k)dk, one 
deals with \{a)da, the elementary intensity due to such point-like sources of angular 
width da. As a result, the visibility is related to the normalized Fourier transform of 
the extended source angular intensity distribution. This problem was already studied 
at the end of the 19th century [O and led Michelson to the construction of his stellar 



interferometer 18 



6. Wigner function in the ordinary channel 

In the previous sections we have proposed a definition of decoherence based on the 
visibility of the quantum interference pattern. As we have seen, this definition has some 
unexpected features, somewhat at variance with expectation. We also found an analogy 
in classical optics. However, alternative definitions of decoherence are possible, based 
on the density matrix and on the Wigner function. Let us therefore briefly recall the 
definition and some properties of the Wigner function. 



The Wigner quasi distribution function |19] can be defined in terms of the density 
matrix p as 

W{x, k) = ^ j di e-'^^x + i/2\p\x - e/2), (51) 



27r 

where x and p = hk are the position and momentum of the particle. One easily checks 
that the Wigner function is normalized to unity and its marginals represent the position 
and momentum distributions 

Trp = J dxdk W{x, k) = 1, (52) 

P{x) = {x\p\x) = J dkW{x,k), (53) 

P{k) = {k\p\k) = [ dx W{x, k). (54) 
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The analyses of the properties of quantum states based on the Wigner function are 
useful because they enable one to make prompt comparisons with fields like quantum 



optics 1 20 1 and quantum tomography |21 



We focus on one-dimensional systems and assume that the wave function is well 
approximated by a Gaussian 

1 



m = im 



(27r52)l/4 

1 



exp 



exp 



[X 



Xo) 



- exp 



(27r52)i/4 

-5\k 



452 



^51 



+ ik^x 
- i{k - ko)xo 



(55) 



k,? 



i{k — kQ)xo 



(56) 



where ip{x) and 0(/c) are the wave functions in the position and momentum 
representation, respectively, 6 is the spatial spread of the wave packet, 6kS = 1/2, 
Xo is the initial average position of the particle and po = hk^ its average momentum. 
The two functions above are both normalized to one. The Wigner function for the state 
(^)- (|56|) is readily calculated 



W{x, k) = — exp 



[x 



Xo) 



252 



exp 



-2s\k-koy 



(57) 



Consider now a neutron wave packet that is split and then recombined in an 
interferometer, with a phase shifter A placed in one of the two routes. The Wigner 
function in the ordinary channel (transmitted component) is readily computed: 



1 



Wo{x, k,A) = — exp[-25\k - kof] 



exp 



(x - Xo + A)' 
252 



+ exp 



+2 exp 



(x - Xo)' 



252 



X 



Xo + 



252 



cos(A;A) 



(5^ 



Notice that, for A 7^ 0, it is not normalized to unity (some neutrons end up in the 
extraordinary channel — reflected component) and that for A = (no phase shifter) one 
recovers (IF/ 



7. Alternative definition of decoherence 



We look at a particular case and assume that the shifts A fluctuate around their average 
Ao according to the Gaussian law (^). The average Wigner function reads 

W{x, k)= J dA w{A) W{x, k, A) (59) 

and represents a partially mixed state. Essentially, this Wigner function represents the 
whole ensemble of neutrons in an experimental run. For the double Gaussian state (|58|), 
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Figure 6. Wigner function in the ordinary channel ( |60[ ) for different values of the 
standard deviation a in (|30|). From left to right, (7 = 0, 0.9, 1.8A. We set xq — 0, 
fco = I.TA^^, 6 — l.lA, Aq = 16. lA. Position x and momentum k are measured in A 
and A~^, respectively. Notice the strong suppression of interference at large values of 
momentum, both in (b) and (c). The interference term in ( |60| ) depends on x and the 
oscillating part of the Wigner function is bent towards the negative x-axis. 



obtained when a neutron beam crosses an interferometer, the average Wigner function 
in the ordinary channel reads 

exp[-2S\k - koY] 



X < exp 



X 



52 



52 + (t2 



exp 



{X + Aq) 
'2(52 + (j2 



+2. 



exp 



(x + 



2(52 



+ 



cos k 



252A0 



xa 



2 52 + 




(60) 



where we set xq = for simphcity. Its momentum marginal (p^ ) (momentum 
distribution function) can be computed analytically and is of interest, because it displays 
fragility at high momenta |TI|, |^: 



P{k) = 1/ — exp \-26^{k 
27r L 



1 + exp 



^2(72 



cos(A;Ao) 



(61) 



The average Wigner function (^) is shown in Figure |^. One clearly observes a 
strong (exponential) suppression of interference at high values of k. Notice that the 
oscillating part of the Wigner function is bent towards the negative x-axis. This is 
due to the x-dependence of the cosine term in (|60D that entails different frequencies for 
different values of x. 

The loss of quantum coherence is clearly visible in Figure ^ as the level of noise 
a increases. One can try to corroborate this qualitative conclusion by introducing a 
quantitative notion of decoherence based on the Wigner function; however, as we shall 
see in a while, one runs into the same kind of difficulties encountered in Sec. ^. We first 
recall that there is an interesting relation between the square of the Wigner function 
and the square of the density matrix: 

Trp2 



dx dk W{x, ky 



2tt 



(62) 
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Figure 7. Decoherence parameter vs coherence length of the wave packet 6 and 
standard deviation of the fluctuation a for a double Gaussian wave packet ( p8| ) in 
the ordinary channel of a MZI. We set fcoAo — 27.4. The decoherence parameter 
is not a monotonic function of a for every value of 6. Notice that e never reaches 
unity (e < 3/4): this is due to the fact that only one Gaussian (in one branch of the 
interferometer) undergoes statistical fluctuations (see Eq. (|60|) and Figured). 



It is therefore possible to define an alternative decoherence parameter [Q, that takes 
into account the coherence properties of the neutron ensemble 

Trp^ 2TTjdx dk W{x,k f 

(Trp)' (^JdxdkW{x,k)) 

This quantity measures the degree of "purity" of a quantum state: it is maximum 
when the state is maximally mixed (Trp^ < Trp) and vanishes when the state is pure 
(Trp^ = Trp): in the former case the fluctuations of A are large and the quantum 
mechanical coherence is completely lost, while in the latter case A does not fluctuate 
and the quantum mechanical coherence is perfectly preserved. The parameter (^) 
was introduced within the framework of the so-called "many Hilbert space" theory of 
quantum measurements and yields a quantitative estimate of decoherence. The 
related quantity Trp — Trp^ was first considered by Watanabe p in 1939 (!). A quantity 
related to e was also introduced in order to get a quantitative estimate of information 



for a quantum system ^2 



It is also worth noticing that the notion of decoherence just introduced is based on 
the square of the density matrix (or Wigner function) and therefore is not accessible to a 
direct measurement procedure. In this sense, it is less "operational" than that discussed 
in Sec. |[ 

The decoherence parameter (|63D is shown in Figure ^ as a function of the coherence 
length of the wave packet 6 in (|55D -(p8D and the standard deviation of the fluctuations 
a. It is not a monotonic function of a for all values of 6. Once again, like in Sec. ^, there 
are situations in which a larger noise yields a more coherent wave packet (according to a 
given definition). The behavior of e has a nontrivial dependence both on the fluctuations 
(cr) and on the wave packet properties (/cq and 6). 
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8. Entropy 

The conclusions of the previous sections can be corroborated and put on a somewhat 
sounder basis by computing the entropy of the distribution of the shifts according to 
the formula 

S= - f dA w{A) logw(A). (64) 



This quantity yields an estimate of the collective "degree of disorder" of the distribution 
of the shifts w{A). One can draw general conclusions about the behavior of S* as a 
function of a parameter a characterizing the width of the distribution. Indeed, let 
w{A;a) be the symmetric distribution with the properties (||), a being its standard 
deviation. By assuming that the distribution function w depends only on the single 
dimensional parameter a, then it must scale according to 

w{A;a) = . (65) 

Therefore 



S{a) = - J dA w{A; a) logw(A; a) 



dA 




log 






— w 




—w 




a 











= - j dA' w{A'- 1) logw(A'; I) + j dA' w{A'- 1) log a 

= S{l) + \oga, (66) 

where 5(1) is independent of a and depends only on the form of the distribution function. 
S{a) is clearly an increasing function of a. 

For example, the Gaussian distribution (BO) yields 



5(a) = loga+llog(27re), (67) 

while the "sine" distribution (|3^) yields 

S(a) = loga-ilog2. (68) 

Therefore, the behavior of the decoherence parameter e as a function of the entropy 
S of the shifts is qualitatively equivalent to its behavior as a function of the standard 
deviation a. Indeed, Figs. ^ || and ^ would differ only for a logarithmic scale on the 
abscissae. As we have seen in this article, in general, the two quantities S and e do 
not necessarily agree: in other words, the loss of quantum mechanical coherence is 
not necessarily larger when the neutron beam interacts with fluctuating shifts of larger 
entropy. 



9. Conclusions 



We have introduced and discussed some interference experiments that display some 
"anomalies" both in the classical and in the quantum domains. The neutron beam 
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partially looses its quantum coherence as a consequence of the fluctuations of the phase 
shifts A. One should emphasize that we have considered the case of "slow" fluctuations, 
in the sense that each neutron crosses a phase shifter of length L, but the length of the 
shifter varies for different neutrons in the beam (different "events"). We have supposed 
that every neutron undergoes a shift A that is statistically distributed according to a 
distribution law w{A). 

We focussed our attention on two alternative decoherence parameters. The first is 
defined in terms of a generalized visibility of the interference pattern in a double-slit 
experiment (MZI) and is more operational. The second hinges upon less operational 
concepts, such as the square of the density matrix. 

All our results corroborate the ideas expressed elsewhere 0] and make it apparent 
that the concept of loss of quantum mechanical coherence deserves clarification and 
additional investigation. It would also be interesting to discuss analogies and differences 
with conceptual experiments in which decoherence is complemented by Welcher-Weg 
information [E^ . 
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